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Abstract 

Generalizing methods developed by Pinn, Pordt and Wieczerkowski 
for the hierarchical model with one component (A^ = 1) and dimensions d 
between 2 and 4 we compute 0(iV)-symmetric fixed points of the hierarchi- 
cal renormalization group equation for some N and d with < d < 4 and 
—2 < N < 20. The spectra of the linearized RG equation at the fixed points 
are calculated and the critical exponents v are extracted from the spectrum 
and compared to Borel-Pade-resummed e-expansion. 
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1 Introduction 



Euclidean quantum field theories and systems of statistical mechanics are 
commonly defined by a generating functional of correlation functions 



where d^i^ is a suitable Gaussian measure with covariance v on the space of 
fields under consideration. The Gaussian measure describes the free theory 
and describes the interaction and will be called the "Boltzmann factor" . 

Let us consider a model with N local degrees of freedom. Famous exam- 
ples of such models in statistical physics are the Ising model with = 1 and 
the Heisenberg model with N = 3. These models correspond to c/)'^ theories 
for fields with real components in the euclidean QFT setup. Here we are 
using a local variable ip G to describe such situations. We analytically 
continue our calculations to noninteger and negative values of A^. The case 

= is an example of physical significance of this continuation, as it is 
used for the description of polymers. 

A^-component models have recently been reconsidered in the literature. 
In this paper we are doing so using the high accuray scheme of the algebraic 
renormalization group. Instead of fine tuning critical parameters it is pos- 
sible to calculate directly the fixed points by using Newtons method. The 
critical exponents can then be calculated from the spectrum of the linearized 
RG equation at the fixed point. We give a list of scenarios for various values 
of N and d. 

We will consider hierarchical models in the sense of I?], Q as an approx- 
imation to "full" field theoretic models models by replacing the covariance 
V = (— A + m?)~^ by a translationally non- invariant hierarchical approxi- 
mation of the massless covariance. They can also be regarded as a system 
of statistical physics with a peculiar structure. The hierarchical covariance 
simulates the long distance behaviour of the "full" covariance and preserves 
the locality of the Bolzmannians under renormalization group transforma- 
tions (RGT). A Boltzmann factor is called local, if it factorizes 



The hierarchical RG equation is a RGT of Boltzmannians Z = e"^ : 





-V(^) ^ TT ^-VMx)) _ 
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For A^-component fields the hierarchical RG transformation is given by 
the non-linear integral transformation 



nz{^) = J d/.^(i_^2)(v)z^ (V + P^) , (2) 

with the Gaussian measure dfiy(ip) = e dipi . . . dip^, ^ > 0, f3 = 

L^'t One may consider this equation as our starting point. A derivation 
of equation (^) for = 1 can be found in 0]. The derivation for general N 
can be done in an exactly analogous way. 

In this paper we present the calculation of the 0(A^)-symmetric fixed 
points and their spectrum. Fixed points of the RGT determine the long dis- 
tance behaviour of the system and are a tool to investigate the continuum 
limit and the thermodynamic limit of such theories. Recently also the stabil- 
ity of 0(A^)-symmetric fixed points of this equation under disturbances with 
cubic symmetry was investigated both for full models and in the hierarchical 
approximation 

The critical exponents are functions of L in the hierarchical approxima- 
tion but they do not vary very rapidly with L unless L is close to 1. The 
limit L ^ 1 is perhaps singular. KoCH and Wittwer have shown the 
existence of a fixed point of this equation with the values N = l,d = 3 and 
L'^ = 2 Q. For L'^ = 2 the model is equivalent to Dysons version of the 
hierarchical model. Also a product of functions can be defined by 



(/ • gm ■■= J d/.^(i_/32)(v)/(V' + P^)gii^ + M • (3) 

The hierarchical RG equation is then given by 

nz = z- z (4) 

and the linearisation L-ji of TZ at the "point" / is given by 

Ln{f)Z = 2f-Z. (5) 
Note that the product defined in this way is commutative but not associative 

It is easy to see that there are two trivial fixed points of the transfor- 
mation TZ (apart from Zq = 0, which does not define a theory), namely 
Zuv = 1 and a Gaussian fixed point 

2/3^-1 2 

ZHTif) = (2/32)e47(i-/3^)^ . 
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Zuy belongs to the massless free theory and Zut to the massive free theory. 
One important feature of the RG equation is 

T^/3,^{ZhtZ) = Zht{T^I3' a''^) 

with f3' = L~^~i = and 7' = L~^7 = 2~'^^. The RG equation with 

/3',7' is called the high temperature picture, because the HT fixed point is 
splitted from the Boltzmannian. In the HT picture Z = 1 belongs to the 
HT fixed point and the inverse of Zht belongs to the UV fixed point. 
The linearisation at Zjjv is given by 

{Ln{Zuv)fm = J dfi^^,_f,20)lf{^p + f3^) . (6) 

Normal ordering is in a sense inverse to Gaussian integration. Using this 
one is able to compute the spectrum of L-jilZuy) analytically. This will be 
done in a first section. After computing the spectrum we go on to discuss 
"algebraic" RG equations and their approximate solution using a computer. 
Another way to compute fixed points is by e-expansion. The e-expansion 
is presumably not convergent. We use Borel-Pade summation of the first 
terms of the series upon the hypothesis that the series is at least Borel 
summable. We conclude with a comparison of our results with those found 
in the literature. 



2 The spectrum at Zuv 

Since the massive theory is not critical the massless fixed point is more 
interesting. The eigenfunctions at this fixed point are - as in the case N = 1 
- Hermite (normal ordered) polynomials. Let Hm be the mth Hermite 
polynomial given by 



^ ' dip"" j^^{m-2j)\23jV ^ ^' 



and 



H^l){^n\ — W^iJ - V (-^)^"^' m-2j 



Note that Hin'if) = <f'^. 
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The A?^-component covariance is diagonal in color space. For a multiindex 
fi = (mi, . . . , mjv) define 7^}7^ by the A?^-fold tensor product of the i?^. 

n^;\^,, ...,iPN) = H^:i{^r) . . . i^it (v^^) = • • • ^Tn ■ 

The normal ordered 0(A'^)-invariant products {(p^)"^ are then given by 



u\+...+VM=n 



UJ^T ^ ' 



;/i+...+i/jv=n 

In the following we will use not use the notation :((^)":-y. Since normal 
ordering is the inverse of Gaussian integration in the following sense 

and 

= mi + . . . + m^r) one finds that the hn^ are a complete set of eigen- 
functions of L-ji{Zijv) in the space L2(M'^, d/x^)*^^^-* of O(A^) symmetric 
functions. The corresponding eigenvalues are An = 2/3^" = 2^~^"'~3r . The 
critical exponent v is calculated by 

_ ln(L) _ ln(2) _ 1 
In(Ai) ~ dln(Ai) ~ 2 ' 

We note that the spectrum in a space of 0{N) symmetric functions is non- 
degenerate and the same as for the case A'^ = 1. 



3 Algebraic Computation 

We now go on to the task of computing non-trivial fixed points and their 
spectrum. For calculations using a computer it is convenient to use the 
so called algebraic RG equation. In this technical section we discuss the 
derivation of this equation. 
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One expands the functions under considerations into a series of Hermite 
polynomials or a power series as 

oo 

or 

oo 
n=0 

These two possibilities to expand the functions minimize the necessary 
work. Consider a general expansion 

oo 
n=0 

where the suitable system of functions. Then 



{nz){^) = 

I d/i^(i„^2)(^)e„(V' + /?¥')em(V' + /3y5) 

n,rTi=0 

oo „ 

n,m,k=0 

For en{ip) = (if^T the expansion of en&m is trivial and on the other hand if 
Cn is a Hermite polynomial the integral is trivial, since Gaussian integration 
is inverse to normal ordering. 

3.1 Expansion in powers of ip^ 

Let us consider first the expansion in powers of (^^. One has 

(7^Z)((^) = Y^ZnZm f dfl.^(^,^p20m + P^fr^^ (9) 
n,m 

= Y.ZnZmht,t~^"^\P^) (10) 



n,m 

„2 



The next step is to express (10) as a sum of powers of (f . To do the necessary 
calculations in a convenient way we introduce the generating function G^'^^ 
of the normal ordered products. 



oo r 



nl 

n=0 
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By reordering of the series for = 1 and induction to general N one finds 

G^-'Ha,^) = L^exp(— ^(^2) 

(1 + 270)- 1 + 270 



Since hn (ip) = (V'^)"' o^is expects G^^\a, tp) = e"'^ which is easily verified. 

As a first step we want to calculate the scalar products {hn\h^})-f in 
L2(M'^, d/^-y). By expanding 

(G(o,-),G(5,.)>7 = / dti^,{^)G{a,^)G{h,^) = {I - A-i^b)-^'^ 



in powers of a and h and comparing coefficients, we find 

n-l 



(/iW, /iW), = n\\{{N + 2i){2j'r5n,m . (11) 



i=0 



Equation (11) can be generalized to compute the coefficients of powers of 
if'^ in equation (|lO|). One calculates 



(G(^')(a, •), GW(6, •))^ = (1 - 2a(7 - 7') - 47^06)-^/^ . 

By expansion into Binomial series two times and comparing coefficients we 
find 

/.(7') .(7)v _ / nto(^^ + 20n!(^)(27T(7-7')'=-" : n<k 

If we want to get the coefficients ch^ in 

00 

4''^ = E-^'^^i'^^ (12) 

n=0 

as in equation (|l^) we can use the last result together with the orthogonality 
of the /i^7^ (111 to find 



ntn(^^ + 2i)C)(7-y)'-" ■■ ri<k 
: n> k 



The derivation is restricted to 7, 7' > 0, but by applying a Gaussian integra- 
tion with covariance 7, the covariances can be changed to 7 — 7 and 7' — 7. 
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We see that equation (|T2|) is correct for all values of 7 and 7'. Inserting the 
coefficients ch^ into (|To|) we get 



00 00 

(n+m) r,2k(,'2\k 



k=Q n,m=0 



By replacing z„ by z„(7(l — ^ that is using the expansion 

(7(1 - 



n=0 



we can eliminate the Term (7(1 — A; £],qj^ ^^j^g final algebraic RG 

equation. 

00 00 

k=0 m,n 

V ' 

(7ez)fc=:4 

rn.+m— 1/A7 , o„-\ 



SriN) :-- 



: k > n + m 



The algebraic RG equation is a mapping of a subset of to instead 
of an integral equation. It is given by 

Z = {zo,Zi,Z2, ...) ^ TZ{z) 

with 

00 

Zk ^ n{z)k = f''Y.^r{N)znZ^ . (13) 

The dimension d is only a numerical parameter in the integral equation. But 
equation ([l^ ) contains N also only in numerical form. We can generalize 
to noninteger values of and investigate other values of N which are of 
interest, e.g. A'^ = 0. This is also possible for the integral equation by inte- 
grating out the 0(A^)-invariant part. The resulting equation is an integral 
equation with only one integration containing a Bessel function and as a 
parameter. 
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3.2 Solving the equation numerically 



To compute approximate solutions using a computer one simply truncates 
the equation. We have to solve a finite system of nonlinear equations of the 
type ( |l3|) . To avoid changing signs of the coefficients in the expansion of the 
nontrivial fixed point we do the calculations in the HT picture. Numerical 
experience suggested as in to use a different normalization, namely 

oo 

The system of equations which was solved numerically was the one consisting 
of Imax + 1 equations (A; € {0, ... , Imax]) given by 

^max 

Zk^Zk-P'^' Sr{N)ZnZm 
n,m=0 

and 

{n + m — 1 / (ouw 

: k > n + m 

This was done by the procedure C05PBF of the NAG library, which is a 
combination of the newton method and the gradient method to find zeros of 
a set of nonlinear equations. The eigenvalues of the linearized RG equation 
have then been computed by the procedure F02AFF. We first reproduced 
the results of [R] to check the correctness of the implementation. 



3.2.1 Effects of the truncation 

For = l,d = 3 an approximation of this type converges for Imax oo 
^ . We investigated the effects of the truncation for general N and there is 
evidence that our approximation converges also. We found it necessary to 
increase Imax with N, see table ||. For further discussion of the effects of the 
truncation see appendix 



4 Numerical results 

We present here the observations made during the numerical investigation 
of the model. There is a bifurcation scenario of the hierarchical RG 0, P], 
which can be understood by means of the implicit function theorem. At 
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^max 


TV = 3 


TV = 10 


TV = 20 


10 


0.77308443210437 






20 


0.76113986112499 


0.92746045611772 




30 


0.76113984902214 


0.91861185755319 


0.98467223300539 


40 


0.76113984902214 


0.91861145661936 


0.96069946588199 


50 




0.91861145661842 


0.96068060650581 


60 






0.96068060512759 


70 






0.96068060512757 



Table 1: Effects of the truncation: as a function of TV and I 



d = 4 a double-well fixed point bifurcates from the UV fixed point, which 
exists for d < A. Generally at d^ = -^zj: k > 2 a A;- well fixed points bifurcates 
from the UV fixed point, which exists for d < d^. This can be understood 
more physically by observing that at dk the eigenvalue Xk of the linearized 
RG equation at Zjjv becomes relevant (Afc > 1). This behaviour has indeed 
been reproduced by our numerical results for k = 2 and for some cases of 
A: = 3. 

For plotting the fixed points one has to choose the parameter 7'. We 
have used the same convention as 

with 7' = ^- The plots show the potentials {V = — ln(Z)) of the fixed points 
in the UV picture normed by 

Z(0) = 1 . 

0(iV)-invariance allowes us to use only one variable for integer TV. We define 
fixed points for noninteger TV by the formula above. 

Close to the bifurcation points our truncation scheme becomes bad due 
to the flatness of the non-trivial fixed point. We present here the conjectured 
behaviour of the system, abstracted from the numerical results. The fixed 
point which bifurcates at d = 4 from the UV fixed point has shown the 
following behaviour. 

• For N = —2 the potential is a single well and the critical exponent is 
0.5. (Fig. 1^, III) The potential could be fitted to an approximatively 
quartic function ip"" (a = 3.982(2)) for 99 = 0.01 ... 0.5 at d = 3. First 
order e-expansion results in v{ip) = ci -1-0293^. 
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For N > —2 the potential is a double-well. The minimum becomes 
deeper with growing N and shifts to bigger arguments. The critical 

exponent i' > ^ grows with N. Zi' {^/N-) converges on a function 
as iV ^ oo. (Fig. I I ) 

u as N ^ CO and ^ ^ as d 4 (Fig. |) . 

N > 1 (Fig. I, D 

— The fixed point exists for 2 < d < 4. 

— u diverges as d \ 2, because the first fixed eigenvalue at the fixed 
point converges on 1: Ai ^ 1 in 



log(Ai)d • 

< iV < 1 (Fig. I, D 

— The fixed point exists for d! < d < A with < < 2. 

— d' grows (non linearly) with from to 2. 



log 2 
log(Ai)d- 



— V diverges as d\ d' , because of Ai ^ 1 in i/ 
. iV = 0. (Fig. I I, 0) 

— The fixed point exists for < d < 4. 

— u diverges as d \ 0, because of ^ — > oo in = j^^^- The first 
eigenvalue converges Ai — > 2. 

• -2<N<0. (Fig. I, 0) 

— The fixed point exists for < d < 4. 

— Ai diverges as d \ 0. There is a < d' < 2, with Ai > Aq = 2 
for < d < d' . d' decreases (non linearly) with growing N from 
2 to 0. 

— The limit of as d \ exists, since the divergencies of Ai and ^ 
cancel. 

• At d = 2 ly diverges as iV / 1, because of Ai 1. (Fig. |, |, 0, 1 ). 

With a view on the theorem of implicit functions Ai ^ 1 indicates pos- 
sible points of further bifurcations or the end of existence of solutions. We 
did not find other solutions below the corresponding values of d and is was 
also not possible to follow the fixed point below these values of d. Based 
on these observations the conclusion is, that these are the points where the 
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existence of the nontrivial fixed points ends. A similar conclusion seems to 
apply for the case d = 2,N 1. The fixed point which exists below = 1 
does not exist for values of A^ > 1. In this case other fixed points are known 
to exist at A^ = 1, d = 2 g. 



A^ = 0.01 


N = 0.1 


A^ = 0.25 


N = 0.5 


N = 0.75 


A^ = 0.9 


0.66(2) 


1.14(2) 


1.46(2) 


1.75(2) 


1.91(2) 


1.97(2) 



Table 2: Values d' for which Xi{d') — 1. 



Furthermore we investigated 3-well fixed points and found the following 
behaviour. 

• There exists a 3-well fixed point for A^ = 1 and 2 < d < 3. (Fig. |8| and 
the figures in Q) 

• There exists a 3-well fixed point for A^ > 1 and d'{N) < d < 3, 
2 < d'{N) < 3. It is unknown whether d' > 2. (Fig. ^ 

• There exists a 3-well fixed point for A^ < 1 and (i'<(i<3,0<d'<2. 
(Fig. I, § 

• The second eigenvalue A2 of the linearized RG equation at the 3-well 
converges on I as d\d'. (Fig. ||) 

5 e-Expansion 

Following closely the lines of one can generalize the e-expansion to the 
model with A^ components {e = 4 — d). We work in the UV picture because 
of the simple structure of the UV fixed point, 

00 

= 1 = J^Z^/iW . 
fc=0 

This results in Zn = (5n,o- Furthermore we expand in terms of hermite 
polynomials since they are the eigenfunctions at ZjjVj see section ||. Again 
the first step is to calculate the algebraic RG transformation for this case. If 
we know the expansion of a product of two hermite polynomials we only have 
to integrate a hermite polynomial, which is trivial. We have to calculate the 
coefficients C^™" in 

00 

k=0 
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This is again done by the method of generating functions. Using the orthog- 
onahty equation (pA|) we get 



fc-i 

i=0 



This leaves us with the calculation of the scalar product. We calculate 

(G(a, ■)G{b, •), G{c, ■))^ = (1 - 47^(06 + bc + ca) - Idj^abc)'^/^ 

and do an a bit lengthy expansion of the right hand side of the equation to 
compare coefficients of o, b, c to find 



(q — nVJq — m)!(g — k)\(n + m + k — IqV. 

n . m , fc < q 

2q<n+m+k 



The final result is 



. ^ — n)!((7 — m)!((7 — A;)!(n + m + A: — 2(jf)! 

2(;<n+m+fe 

It was verified using Maple that for A'^ = 1 these coefficients are the same as 
those found by |Q] . The integration of h^^J^ in the RG equation is now trivial 
and only amounts to a multiplication with /J^*^. 

To make the algebraic RG equation again independent of 7 we use the 
expansion 

00 

z = 5]z„7-"/i^7^. 

n=0 

Inserting this expansion into the RG equation results in 

00 00 

k=Q m,n=0 

^ V ' 

(7^Z)fc=:4 

We define the 7-independent structure coefficients C^"^ by 

/-ynm , pnm^—n—m+k 
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and get the following algebraic RG equation z = (zq, zi, . . .) i-^ TZz with 

oo 

Zk^{nzh=P^^ Cr^nZm- (15) 

n,m=0 

5.1 e-expansion of the fixed points 

The e-expansion is an expansion of the coefficients zi of a fixed point into 
an asymptotic series of the form 

oo 

We write the fixed point equation in the form 

n,m 

Expanding the left and the right hand side of the equation in powers of e 

(k) 

yields recursive equations for the z^ . The calculation of these equations 
does not differ in an essential way from the one in and we will leave it 
out here. One only has to replace the structure coefficients by the C^™ we 
just calculated. 

To first order we find for the fixed point 

In 2 



16(iV + 8) ' 2 



(7) 

For we have 



We note that for N = —2 and e > the fixed point has the form of a quartic 
single- well and for > —2 it is a double- well. 

Figure ^ presents some fixed points calculated to order 1, 2 and 3 by 
e-expansion compared to numerically calculated fixed points. 

5.2 e-expansion of eigenvalues 

Given the fixed point we can calculate the linearisation of the RG equation 
at this fixed point. The linearisation is given by the matrix A(z) 



mn 



m=0 
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The eigenvalue equation is then 

A(z)v = Xv . 

We get the e-expansion of A{z) by inserting the expansions of /3^' (with 
coefficients b^^^) and Zm 

oo k 

A(z),. = 2'--^Y.ET.crmz^M'r^e>' ■ m 

k=0 u=0 m 

Expanding also the eigenvectors 



k=0 

and the corresponding eigenvalues 



k=0 

we can use a generalisation of nondegenerate perturbation theory. For details 
we refer to ||^, ^. Summarizing we get the following equations (/=„ = 2) 

A{zfj = 2'-^Y.E criN)z^:^b^':r^ (17) 

vj^^ = {Sn,^^) and = 2^'^ , ^ € No (18) 

oo 

= (^?),<')) = E'^^."-S = *i (19) 

n=0 

k—l oo oo k—l oo 

>-p = E E e E a:;,-"'";',:; = e e as-"'"!.:? . <») 

n=Om=0 i=0 n=0 j=0 

^ fc— ^1 / oo \ 

« = E E (a?^-'^)-^) - Af-^.S (21) 

2(2 -2-T7)^^\^^y M, ; M MPy 

With these equations it is possible to calculate the e-expansion of the eigen- 
values and eigenvectors recursively order by order. To each order the sums 
are finite. 
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The critical exponent v belongs to Ai. The first order coefficient is 



In 2 / 1 + 2 



Summing up to first order 



V2{l + 



In 2 / 1 N + 2 



N + 8 



and expanding the logarithm in ^ = in 2"'''^ S®^ 



2N + 



so that 



1 lA^ + 2 ^,2, 

2 + 4ivT^^ + ^^^)- 



The e-expansion was calculated up to 6th order using Maple V Release 3. 
The results of Q for the case = 1 were reproduced. In appendix some 
analytical results of the e-expansion are given. 



d 


k = 1 


fc = 2 


it = 3 


fc = 4 


fc = 5 


numerical 


BP 


3.9 


0.5065 


0.5066 


0.5066 


0.5066 


0.5066 


0.50659 


0.50659 


3.8 


0.5135 


0.5140 


0.5137 


0.5139 


0.5137 


0.51380 


0.51380 


3.7 


0.5211 


0.5222 


0.5211 


0.5221 


0.5208 


0.52157 


0.52157 


3.6 


0.5293 


0.5313 


0.5286 


0.5319 


0.5261 


0.52990 


0.52988 


3.5 


0.5382 


0.5412 


0.5359 


0.5442 


0.5266 


0.53877 


0.53873 


3.4 


0.5477 


0.5521 


0.5430 


0.5602 


0.5176 


0.54820 


0.54811 


3.3 


0.5580 


0.5641 


0.5496 


0.5821 


0.4945 


0.55821 


0.55761 


3.2 


0.5691 


0.5771 


0.5555 


0.6128 


0.4557 


0.56884 


0.56841 


3.1 


0.5810 


0.5913 


0.5607 


0.6570 


0.4045 


0.58012 


0.57946 


3.0 


0.5940 


0.6068 


0.5651 


0.7236 


0.3488 


0.59209 


0.59040 



Table 3: i'{N = 0) calculated by different methods explained in the text. 



naive 

_ In 2 



din A 



way, that is by summing the e- 
■. This is called "naive", because 



We calculated 1^ first in the 
expansion of A and calculating i' 
the e-expansion is presumably not convergent. Upon the hypothesis that it 
is at least Borel-summable, we calculated the Borel sum of the expansion 



max. order > (^]^^ 

k=0 
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Then we calculated the diagonal Pade approximation Q(e) of this sum and 
did a numerical Borel transform of Q 

/•oo 

Ai(e) = / Q{te)exp{-t)dt . 
Jo 

Table ^ shows values of the critical exponent calculated by "naive" 
summation of the e-expansion up to order k, by numerical computation of 
the eigenvalues and by Borel-Pade summation of the e-expansion up to 6th 
order. 

Table ^ shows values of u calculated by Borel-Pade approximation of 
the e-expansion up to 6. order. For small e there is a very good agreement 
between the values calculated by solving the truncated system and those 
calculated by Borel-Pade summation of the e-expansion. If we go to d = 3 
only 2 to 3 digits agree. Unfortunately we do not have error estimates for 
the Borel-Pade summation. 



6 Comparison with other results and Conclusion 





TV = -1 N ^0 


N = 1 


TV = 2 


TV = 3 


TV = 4 




0.541536 0.592086 


0.649570 


0.708225 


0.761140 


0.804364 










0.64957 












Pol 




0.6496 


0.7082 


0.7611 


0.8043 






WH 


0.6066 


0.6895 


0.7678 


0.8259 


0.8648 




1 


WH 




0.6896 


0.767 


0.826 


0.865 




10 HT 




0.6301(18) 


0.6734(28) 


0.7131(40) 


0.7361(68) 




5| pert. 


0.5882(11) 


0.6304(13) 


0.6703(13) 


0.7073(30) 


0.7410(60) 



Table 4: Comparison of values of v for d = 3. 



Table 3 of [^j contains values of the critical exponent u calculated from 
the Wegner- Houghton and the Polchinski RG equations. It is remarkable 
that our values coincide to three significant digits with those found by the 
Polchinski RG equation. See table ^. The results of the Wegner-Houghton 
equation differ from our results, but are compatible with those in Com- 
pared to results for the "full" model the hierarchical model overestimates 
the critical exponent. We refer to the results of found by convergent 
high temperature expansion and of calculated by perturbation theory. 
For N ^ oo and —2 the hierarchical model has the same exponents 

as the "full" model. 
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We have presented the calculation of 0{N) symmetric fixed points of the 
hierarchical RG transformation by e-expansion and by numerical methods 
based on the expansion of functions. The results of the e-expansion are com- 
patible with those found by numerical computation when e is not to large. 
Furthermore the results are compatible with results found by investigating 
the Polchinski equation. 

Considering this, 0{N) symmetric solutions of the hierarchical approxi- 
mation seem to be quite well understood. Mathematically there is still room 
for improvement since rigorous proofs of the existence of the fixed points cal- 
culated in this paper are to my knowlege still lacking for most values of N 
and d with the exception N = 1, d = 3, = 2. Cassandro and Mitter 
have found other fixed points for d < 2 and > 1 Q. It should be possible 
to find these fixed points by algebraic computation. In this work no fixed 
points for these values of N and d have been calculated. It would be inter- 
esting to know if at least for 2 < d < 4 the bifurcation scenario leads to all 
fixed points of the RG equation. 

Other - more physical - routes from this point are systematic corrections 
to the hierarchical approximation and the study of fixed points with less 
symmetry, e.g. |Q. 
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A Effects of the truncation 

Table ^ shows the first 10 coefficients of the fixed point for N = 3,d = 3. If 
we go to Imax = 40 the first 10 coefficients are the same as for Imax = 30. For 
l^^^ = 40 we found zIq = 3.2209843597592 • lO'^o, z^^ = 1-0554235840436 • 
10^13 and 4o = 5.0374367851155 -10^^ So we can use Imax = 30 for iV = 3. 

= 20 was the greatest value of A'^ used during numerical calculations. 
See table ^ for a selection of coefficients of the fixed point. One can observe 
the increasing difficulty of finding the fixed point because the first 10 coeffi- 
cients have the same order of magnitude and they are quite small. One has 
to go to Imax = 60 before the ffist 10 coefficients are constant under further 
increasing of Imax- Note further that the coefficients of the fixed point do 
not decrease as fast as for lower values of N. 
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2fc 


imax — 10 


imax — 20 


imax — 30 




3.8716018937757 ■ 10"^ 


3.9015294442769 ■ 10~^ 


3.9016294689101 ■ 10"^ 




2.7472709654097 ■ 10"^ 


2.7468216314596 ■ 10~^ 


2.7468216310166 • lO"'^ 




2.0416445882152 ■ 10"^ 


2.0327597434157 ■ 10~^ 


2.0327697360790 • 10"^ 




1.4092248806621 ■ 10"^ 


1.4013502479152 ■ 10~^ 


1.4013502414472 • 10"^ 


^2 


8.9625774854293 ■ 10"^ 


8.9257511848321 - 10^^ 


8.9257511547236 • 10"^ 




5.2733434885004 ■ 10"^ 


5.2748796970937 ■ 10^^ 


5.2748796984464 • 10"^ 




2.8856392240700 ■ 10"^ 


2.9099303116770 ■ 10^^ 


2.9099303316272 • 10^^ 


^? 


1.4747091527680 ■ 10"^ 


1.5072774239603 ■ 10^^ 


1.5072774508427 • 10^^ 


^1 


7.0528086140537 ■ 10"^ 


7.3689239174505 ■ 10"-^ 


7,3689241862291 • 10^'^ 




3.1541868355104 ■ lO"" 


3.4157386267627 ■ 10"^ 


3.4157388689480 • 10~^ 




1.3139603409337 ■ 10~^ 


1.5071112742866 ■ 10"" 


1.6071114902066 ■ 10~^ 














6.0275646988651 ■ 10~" 


5.0374367836766 -10"" 










^31) 






1.0553957852574 ■ 10~^^ 



Table 5: Effects of the truncation for A/" = 3 





Imax = 40 


Imax = 60 


Imax = 60 




6.7135699167269 ■ 10"* 


6.7137831100317 ■ 10"* 


6.7137831236764 ■ 10"* 




5.4722582473628 ■ 10~* 


5.4724091850071 ■ 10~* 


5.4724091946673 • 10~* 




.').22()24122'i4'i2() ■ 10^'^ 


.").22()36.116'i2976 - 10"'' 


".22f):!6.'')17:!2299 ■ 10"* 




r).0324322.:j(i.:jS()9 ■ IIJ ' 


r).():i2.r)33944().:j37 - iO ' 


r),i):!25:i:!!).';iJ56i« ■ iO "* 




4.7793545596989 ■ 10"* 


4.7794355811177 ■ 10~* 


4.7794366863032 10"* 




4.4366461453368 ■ 10"* 


4.4367080895417 ■ 10"* 


4.4367080936062 ■ 10"* 




4.0145889929024 ■ 10"* 


4.0146340398981 ■ 10"* 


4.0146340427811 • 10~* 




3.5381684226392 ■ 10"* 


3.5381992679245 ■ 10"* 


3.5381992698987 • 10"* 




3.0372937145831 ■ 10"* 


3.0373132889402 ■ 10"* 


3.0373132901930 • 10"* 




2.5407316573042 - 10"* 


2.5407428246010 ■ 10"* 


2.5407428253157 • 10"* 


29 


2.0723798082894 ■ 10~* 


2.0723851264509 ■ 10~* 


2.0723861267913 • 10~* 












1.1315164646491 ■ 10"" 


1.1338793621549 • 10"" 


1.1338796346646 ■ 10"" 


^37 


5.7800321411117 • 10~^" 


6.7989994920249 ■ 10~^" 


5.7990021969982 ■ 10"*" 


^38 


2.9151584691444 ■ 10~^" 


2.9299955707748 ■ 10"^" 


5.7990021969982 ■ 10"*" 


2^.3!) 


1.4515862129964 ■ 10~^" 


1.4628783782344 -10"^" 


1.4628809676271 ■ 10~^" 


"In 


7.1353999581577 ■ 10"^^ 


7.2189285198076 ■ 10"^^ 


7.2189630200667 ■ 10"*^ 














3.3911297113494 ■ 10~^* 


3.3943990779867 ■ 10~** 










2fin 






6.0217088064527 • 10~'° 



Table 6: Effects of the truncation for N = 20 
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By increasing Imax to 70 the first 10 coefficients are identical to those 
for Imax = 60. We found for l^ax = 70: z^q = 6.0220916546711 • 10"^^^ 
z^o = 3.3943995424453 • lO'^"^ and z|o = 7.2189530204572 • 10"". 

B Numerical results: the exponent u 



d 


iV = -1 


iV = 


d 


iV = -1 


iV = 









2 


0.572496 


0.766551 


0.1 






2.1 


0.570806 


0.742861 


0.2 






2.2 


0.568709 


0.720999 


0.3 




3.357170 


2.3 


0.566240 


0.700752 


0.4 


0.528092 


2.557642 


2.4 


0.563437 


0.681937 


0.5 


0.534783 


2.094339 


2.5 


0.560333 


0.664404 


0.6 


0.541065 


1.794838 


2.6 


0.556966 


0.648025 


0.7 


0.546882 


1.585792 


2.7 


0.553369 


0.632691 


0.8 


0.552197 


1.431311 


2.8 


0.549577 


0.618308 


0.9 


0.556982 


1.312021 


2.9 


0.545622 


0.604796 


1 


0.561217 


1.216665 


3 


0.541536 


0.592086 


1.1 


0.564886 


1.138306 


3.1 


0.537351 


0.580118 


1.2 


0.567979 


1.066439 


3.2 


0.533094 


0.568842 


1.3 


0.570489 


1.016113 


3.3 


0.528794 


0.558214 


1.4 


0.572419 


0.967166 


3.4 


0.524479 


0.548199 


1.5 


0.573773 


0.924112 


3.5 


0.520177 


0.538766 


1.6 


0.574564 


0.885845 


3.6 


0.515917 


0.529896 


1.7 


0.574809 


0.851529 


3.7 


0.511728 


0.521574 


1.8 


0.574528 


0.820526 


3.8 


0.507645 


0.513801 


1.9 


0.573748 


0.792332 


3.9 




0.506593 



Table 7: Numerically calculated values of v 



Table |^ and table ^ show some numerically calculated results of the 
critical exponent v for different values of N and d. 

C Some results of the e expansion 

In this appendix we present some of the results of the e expansion calculated 
using Maple 

C.l e expansion of the coefficients Zm up to 2. order 

T = V2 und R = ln2 

4°' = 1 
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d 


N = 1 


N = 2 


N = 3 


N = 5 


iV = 10 


2 












2.1 


2.088862 


9.193920 








2.2 


1.362344 


4.045225 








2.3 


1.099159 


2.288147 


2.889098 


3.126743 




2.4 


0.957038 


1.559028 


2.030739 


2.287148 




2.5 


0.865337 


1.214032 


1.535818 


1.782761 


1.908031 


2.6 


0.799849 


1.022244 


1.236507 


1.449465 


l..")74435 


2.7 


0.719927 


().l)0122() 


1.0i()."):'>0 


1.2i7()8:i 


i.:>:57()54 


2.8 


0.710110 


0.817609 


0.918937 


1.051532 


1.160436 


2.9 


0.677294 


0.755940 


0.828437 


0.929516 


1.024914 


3 


0.649570 


0.708225 


0.761140 


0.837755 


0.918611 


3.1 


0.625700 


0.669949 


0.709101 


0.767069 


0.833820 


3.2 


0.604841 


0.638382 


0.667563 


0.711338 


07652475 


3.3 


0.586400 


0.611780 


0.633545 


0.666445 


0.699111 


3.4 


0.569946 


0.588980 


0.605103 


0.629584 


0.662617 


3.5 


0.555163 


0.569176 


0.580926 


0.598813 


0.623711 


3.6 


0.541815 


0.551801 


0.560104 


0.572765 




3.7 


0.529725 


0.536443 




0.546388 




3.8 


0.518769 


0.522810 


0.526130 


0.531208 




3.9 


0.508867 


0.510699 


0.512199 


0.514498 





Table 8: Numerically calculated values of u 







(2) 


(0) 

1 

(1) 

1 

(2) 



(0) 
2 

(1) 



(2) 



(0) 
3 

(1) 

(2) 
3 

(0) 



-NR^{N + 2) 
32{N + 8)2 




8(r- 2)(iV + 8)2 


R 

16(Af + 8) 

r( - 8{N + 8f + 4:{-N^ + UN + 200)R + 6{N + SfT + 



{-92N - 424 + 3N'^)TR^ / (^{N + 8)^(-384r + 512) 



R^ 



32(Af + 8)2(T- 1) 



= 



20 



(2) 



= 



512(iV + 8)2 



C.2 e expansion of the exponent u 



(1) ^ l N + 2 
4iV + 8 

16(7V + 8)^(4 -ST) 
(^~8iN + 8){N + 2) + 12(7Af + 20)R 

6{N + 8)(iV + 2)T - 7{7N + 20)RTj 

(3) ^ (iV + 2) 

32(7V + 8)5(4 -3r)2 

^ - 68(iV + 8f{N + 2f + 45(Af + 8){7N + 20)(5iV + 16)i? 

(681iV^ - 440Af^ - 2648Af + 42016)i?^ + 48(iV + 8f{N + 2.fT 

- 32(Ar + 8){1N + 20)(5iV + 16)i?r 

- 6(79iV^ - 60Af^ - 432Af + 4544)i?^T ) 



The zeroth and first order of the e expansion are the same for the hier- 
archical and the full model. For the full model we have 



(2) _ N + 2 ,^^2 

+ 

~ 32{N + 8)5 



and 

lyf^ = ;j^^^^-^(2iV^ + 89iV3 + 1412iV2^ 5904iV 
8640 - 192(5iV + 22) (iV + 8)t) 



with t ^ 0.60103 [T^. Note that the prefactors (^^^^yi and (^^^^y^ are 
the same in both expansions and that the limit ^ oo is the same for 
both models and equals the coefficients of the expansion of = ^j^^- 
(^«^i,..(2)^iandr.(3)^^) 
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JV 




3.8 




3.7 




3.6 




3.5 




3.4 




3.3 




3.2 




3.1 




3.0 


— 2 





50000 





50000 





50000 





50000 





50000 





50000 





50000 





49996 





49994 


— 1 





50764 





51172 





51586 





51995 





52371 





52744 





53299 





53266 






() 





51380 





52157 





52988 





53873 





54811 





55761 





56841 





57946 





59040 


1 





51877 





52973 





54183 





55525 





56972 





58789 





60133 





62597 






2 





52281 





53644 





55181 





56920 





58817 





61332 





63773 





67734 






3 





52613 





54199 





56010 





58093 





60520 





63391 














4 





52889 





,54660 





56699 





59064 





61836 





65126 





69084 





73922 





79939 


5 





53121 





55046 





57274 





59872 





62928 





66562 





70931 





76246 





82801 


6 





53318 





55372 





57757 





60546 





63834 





67747 





72449 





78161 





85185 


7 





63487 





55650 





58167 





61113 





64589 





68727 





73696 





79724 





87121 


8 





53633 





55890 





58517 





61593 





65224 





69543 





74724 





81000 





88690 


9 





53761 





56098 





58818 





62004 





65761 





70227 





75578 





82050 





89964 


10 





53873 





56279 





59079 





62357 





66220 





70806 





76293 





82919 





91008 


11 





53973 





56439 





59308 





62664 





66615 





71300 





76897 





83646 





91868 


12 





54061 





56581 





59510 





62933 





66958 





71725 





77412 





84258 





92585 


13 





54141 





56707 





59688 





63169 





67258 





72093 





77854 





84778 





93187 


14 





54212 





56821 





59847 





63378 





67522 





72415 





78237 





85225 





93697 


15 





54277 





56923 




59990 





63565 





67755 





72698 





78572 





85611 





94133 


16 





54336 





57015 




60119 





63732 





67963 





72949 





78866 





85948 





94510 


17 





54390 





57099 




60235 





63883 





68150 





73172 





79126 





86243 





94838 


18 





54440 





57176 




60341 





64019 





68318 





73372 





79357 





86504 





95125 


19 





54485 





57246 




60438 





64143 





68469 





73551 





79564 





86737 





95379 


20 





54527 





57311 




60526 





64256 





68607 





73714 





79751 





86945 





95605 


1 

d-2 





55556 





58823 


0.62500 


0.66667 


0.71429 


0.76923 


0.83333 


0.90909 


1.00000 



Table 9: v: Diagonal Borel-Pade summation of the e-expansion up to 6. 
order. 



D Figures 
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Figure 1: Potentials of the fixed points {N = 0,d = 0.3 to 3.9, steps of jq). Deeper 
potentials belong to lower values of d. 




Figure 2: Potentials of the 3-well fixed points {N = 0,d= 2.0 to d = 2.9). The 
local extremum grows with decreasing d. 
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OS a function of d and N 



Figure 3: The critical exponent v for < < 4 and (from the lower left to the 
upper right) N = -1, -0.5, 0, j, \, |, ^, 1 to 5, 10 and 15. The plotted line 
is the asymtotic value oiv iov N ^ oo 



V [d = 2, N = -2 to 1) 



Figure 4: The critical exponent v ioT: d = 2 and = — 2 to 1. 
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potentiol of the fixed point (d = 2 N = -2 to 9) 



2 



-1 





'1 1 


1 ' 








N = -2.0 and -1.9 




^ 




1 0.0010 






W// 




- 0.0005 










0.0000 










- -0.0005 










" -0.0010 










; 0.000. 100. 200. 300. 400. 5l ''-^M: 


Ipip 













































12 3 4 5 



Figure 5: Potentials of the fixed points {d = 2,N = —2 to 0.9, steps of j^). Deeper 
potentials belong to bigger values of N. 



V„(x) (d = 3 N=1 to 20) 



> 




-0.5 I I I I 

1 2 3 4 5 



Figure 6: Potentials of the fixed points (d = 3, A'' = 1 to 20) after rescaling. The 
potential deviating most from the "limit potential" belongs to A'' = 1. 
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Figure 7: Eigenvalues for iV < 0. The unbroken lines belong to the eigenvalues at 
the UV fixed point, continued to d < 2. The (• • •) lines belong to the eigenvalues at 
the 2-well fixed point. The picture for A'' = —1.5 shows additionally the eigenvalues 

at the HT fixed point ( ) and the picture for = the eigenvalues of the 3- well 

fixed point {■■■), which bifurcates at d = 3 from the UV fixed point. At every fixed 
point we also have a "volume eigenvalue" Aq = 2. 
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eigenvalues at the UV fi-<ed point ond 2 well (N=0.01 ) 




eiqenvolues ot the UV fixed point ond 2 well (N=0 5) 




eigenvalues at the UV fii<ed point ond 2 well(N=0.1) 




eiqenvolues ot the UV fixed 


point, 2- ond 3 well(N = 1) 












Figure 8: Eigenvalues for > 0. The unbroken lines belong to the eigenvalues 
at the UV fixed point. The (• • •) lines belong to the eigenvalues at the 2-well fixed 
point. The picture for A'' = 10 shows additionally the eigenvalues at the N ^ oo 

fixed point ( ) and the picture for = 1 and N = 5 the eigenvalues of the 3- well 

fixed point {■•■), which bifurcates at d = 3 from the UV fixed point. At every fixed 
point we also have a "volume eigenvalue" Aq = 2. 
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Figure 9: e-expansion to 1. ( ), 2. (— • — •) and 3. order (• • •) compared to 

numerically calculated fixed points. 
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